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Turbulent combustion can not be simulated adequately by conventional moment closure turbulence models. 
The difficulty lies in the fact that the reaction rate is in general an exponential function of the temperature, and 
the higher order correlations in the conventional moment closure models of the chemical source term can not be 
neglected, making the applications of such models impractical. The probability density function (pdf) method 
offers an attractive alternative: in a pdf model, the chemical source terms are closed and do not require additional 
models. 

The partial differential equation for the probability density function, P , can be written as 
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where the terms represent the time derivative, mean convection, chemical reaction, turbulent convection, and 
molecular mixing, respectively. The fact that the pdf equation has a very large dimensionality renders finite 
difference schemes extremely demanding on computer memory and CPU time and thus impractical, if not entir y 
impossible. A logical alternative is the Monte Carlo scheme, wherein the number of computer operations increases 
only linearly with the increase of number of independent variables, as compared to the exponential increase in a 

conventional finite difference scheme. 

A grid dependent Monte Carlo scheme Mowing that of J.Y. Chen and W. Kolimann has been studied in the 
present work. In dealing with the convection and diffusion of the pdf, the pdf equation is discretized on a given 

grid, e.g., _ - 

Px+4xj = + A?\P rj + Pxj — l 

where 

a/ + Pi + 7 j = 1 

However, if this is the .only restriction satisfied by the numerical algorithm, the mass fractions may not be con- 
served due to re-contamination, and directional-bias also appears. These phenomena are lUustrated in Figure . 
Consider a mixing layer; use white balls to represent contaminants in the upper stream and black balls to re P r f* en 
contaminants in the lower stream. As the two streams move toward right, the location of the white balls and black 
balls are interchanged randomly to simulate convection and diffusion. From Figure 1, it is dear that directional- 
bias caused by recontamination caused the center of the mixing layer to drift downward. ( Directional-bias can 
be partially corrected by changing sweeping directions. ) One also notices that after the first marching step, the 
conservation law is violated, reflected in the Figure as missing white or black balls. 

It is found that in order to conserve the mass fractions absolutely, one needs to add further restriction to the 
scheme, namdy 

Q] + 7i = Q i-i + T>+i 

A new algorithm was devised that satisfies this restriction in the case of pure diffusion or uniform flow problems. 
Using the same example, it is shown that absolute conservation can be achieved. This result is shown in Figure . 
One can see that" the diffusion process is symmetric, and the problem of directional-bias is dimmaled. 

Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the 
error considerably. 
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Figure 1. Illustration of directional-bias and recontani nation 
problems in Monte Carlo simulation. 
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Figure 2. Solution obtained by using additional constraint 
in the solution orocedure. 


